function np = normalpdf(x,mu,sigma)
sigma = sqrt(sigma^2);
np = 1/(sigma*sqrt(2*pi))*exp(-((x-mu).^2)/(2*sigma^2));
